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Zastosowanie wirującej elektrody dyskowej oraz wirującej elektrody 
dyskowej z pierścieniem do badania kinetyki i mechanizmu reakcji 
elektrodowych. Podstawy teoretyczne 


The Application of the Rotating Disc and Rotating Ring-Disc Electrodes for 
Investigations of the Kinetics and of Mechanism of Electrode Reactions. 
Theoretical Basis 


1. WSTĘP 


Wirująca elektroda dyskowa (RDE) i wirująca elektroda dysk-pierścień 
(RRDE) znajdują duże zastosowanie w elektrochemii. Bierze się to stąd, iż ma- 
ją wiele zalet wyróżniających te właśnie elektrody spośród innych elektrod sta- 
łych. Opracowano teorie analityczne, które opisują podstawowe parametry ele- 
ktromechaniczne i hydrodynamiczne procesów przebiegających na tych ele- 
ktrodach. Stosowanie powyższych elektrod umożliwia prowadzenie pomiarów 
w warunkach stacjonarnych, co pozwala uniknąć wielu rodzajów zakłóceń, ta- 
kich jak niestałość mierzonych prądów, mała powtarzalność pomiarów itp. Nic 
więc dziwnego, że od momentu wprowadzenia tych elektrod do badań elektro- 
chemicznych pojawiło się w stosunkowo krótkim czasie w literaturze fachowej 
wiele artykułów przeglądowych [1-9] i monografii [10—12] dotyczących teorii 
i zastosowania RDE i RRDE. Wynika stąd duże znaczenie tych elektrod dla ba- 
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dań nie tylko w elektrochemii, ale także w hydrodynamice i badaniach zjawisk 
transportu masy. 


2. PODSTAWY TEORETYCZNE METODY RDE 


Wirująca elektroda dyskowa ma postać cylindra z materiału izolacyjnego (np. te- 
flonu), w którego dolnej części znajduje się osadzony centrycznie dysk, stanowiący 
właściwą elektrodę. Materiałem konstrukcyjnym elektrody są zazwyczaj metale (np. 
platyna, srebro, złoto, miedź) lub amalgamaty (Au-Hg, Cu-Hg), a także różne postaci 
węgla (grafit pirolityczny, grafit szklisty, pasty grafitowe). Elektroda dyskowa wiruje 
w roztworze ze Ściśle określoną szybkością, a jej ruch obrotowy wywołuje przepływ 
cieczy. Wynika stąd, że zjawiskami mającymi decydujący wpływ na właściwości 
RDE są zjawiska przepływu cieczy i transportu masy. 

Podstawy teoretyczne metody stworzył Levich [13]. Zostały one ekspery- 
mentalnie zweryfikowane przez Sivera i Kabanova [14]. Podstawowym 
równaniem opisującym transport masy jest równanie" (1) 


dc; /dt = D; * Ac, v grad c, —u,* „v(C,* E), (1) 


gdzie indeks i oznacza i-tą substancję, D; jest jej współczynnikiem dyfuzji, c; stę- 
żeniem, v wektorem prędkości przepływu cieczy, u; ruchliwością i-tej substancji, 
a E oznacza natężenie pola elektrycznego w roztworze. W przypadku przepływu 
laminarnego, jeSli znane są składowe wektora prędkości przepływu cieczy v i przy 
założeniu, że transport migracyjny jest pomijalnie mały (E=0, czyli transport czą- 
steczek obojętnych lub nadmiar elektrolitu podstawowego) powyższe równanie 
można zapisać w następującej postaci: 


dc/dt + v,* dc/dr + (v, /r) + óc/0Q + v,* ðc/ðy = 
=D [@c/ay? + 3? c/8, + (1/r) - dc/ar + (1/r2) - e/g’), (2) 


gdzie v„ v, i v, są składowymi wektora v. Ze względu na symetrię osiową można 
opuścić w równaniu (2) pochodne względem kąta Q: 


dc/dt = D - (9? c/ar? + (1/r) - dc/dr + 8*c/dy*) — v, dc/dr — v, dc/ay. (3) 


* W dalszej części artykułu autor, powołując się na poszczególne równania zapisane w tym 
tekście, oznacza je cyframi arabskimi w nawiasie okrągłym — przyp. red. 
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Składowe v, i v, wektora v można wyrazić następującymi zależnościami [15, 
16]: 


v,=r-@- [0,51 -y—0,51 - y? + 0,205 *y*—...] dla 0< y< d,, (4) 
v, = -0,886 : (wv)'? dlay + ©, (5) 
v, = (wv)!? - [0,51: y? — 0,338 - y? + 0,104 * y*...] dla 0< y< ô, (6) 


gdzie w jest szybkością wirowania elektrody w radianach na sekundę, y = (w/v)? 
Y, © jest grubością warstwy granicznej (warstwa Prandtla), a v oznacza lepkość 
kinematyczną roztworu. Według Levicha [17] | 


5, =3,6: (v/w)'? (7) 


Wartość 6, jest więc niezależna od zmiennej r, tzn. nie zmienia się na powierz- 
chni dysku i maleje wraz ze wzrostem szybkości wirowania elektrody w. Tak więc 
w dużej odległości od powierzchni elektrody strumień cieczy jest skierowany pro- 
stopadle do jej płaszczyzny. W pobliżu jej powierzchni ten przepływ zmienia się 
na odśrodkowy, wskutek czego ciecz dopływająca do środka dysku jest odprowa- 
dzana w stronę jego brzegów. Bezpośrednio przy powierzchni dysku wytwarza się 
warstewka cieczy (zwana warstwą Prandtla), która porusza się wraz z elektrodą. 
W przypadku dysku składowa v, ma decydujące znaczenie, ponieważ substancja 
elektroaktywna przemieszcza się do elektrody tylko w kierunku y. Równanie (3) 
można zapisać w następujący sposób: 


v, dc/dy = È c/dy’, (8) 
zakładając spełnienie następujących warunków: 
1) dc/St=0 (tzn. proces jest stacjonarny); 
2) D jest niezależne od stężenia; 
3) c jest niezależne od r; 
4) powierzchnia dysku jest gładka i jednorodna. 


Równanie (8) można rozwiązać dla odpowiednich warunków brzegowych, tj. 


c=zc? dlay >», (9) 


ip = D-(ac/8y)„, =k-c(0) dla y=0, (10) 
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gdzie j, oznacza przepływ substancji elektroaktywnej na dysku, c° jest stęże- 
niem substancji elektroaktywnej w głębi roztworu, a k — stałą szybkości reakcji 
elektrodowej (pierwszego rzędu), która zależna jest od potencjału elektrody: Po 
odpowiednich przekształceniach [5, 10] otrzymamy następującą zależność wyra- 
żającą wielkość tego przepływu: 


jp = k*c?/l+A,k:C 8 D2) (11) 
gdzie A, =1,288, a stała C = 0,51: w 3? -v-'” . Stąd wyrażenie dotyczące prądu dys- 
ku ma następującą postać: 

i Sèn FA. j= n- EF A-k-ciit+A,k- CD), (12) 
gdzie A oznacza powierzchnię dysku. Przy zmianie potencjału wartość stałej k 


zmienia się od małych wartości (A, kC 13 D-23 << 1) — i wtedy 


= tn Fs Aak co — (13) 


i, = 
do bardzo duzych wartości, gdy A, KC“? D23 >>1. Prąd dysku osiąga wtedy war- 
tość graniczną opisaną następującym wyrażeniem: 

inp =n*A-A-D -c° (A4 DC) = 
=n-F-A-D-c°fS (14) 
= 0,62 -n-F-A- D™ -y- 6. 912. (0, 


gdzie 5 oznacza grubość granicznej warstwy dyfuzyjnej, przy czym: 


S=4,-D*B = 1,61 - (DA): (v/w)'2, (15) 


Zmiana postaci równania opisującego natężenie prądu dysku w zależności od 
wartości k (wzrastającej wraz ze wzrostem potencjału) wyraża przejście od prądu 
dysku prądu kontrolowanego kinetyką reakcji elektrodowej (13) do prądu kontro- 
lowanego przez transport (14). Parametr 5 dzieli przestrzeń pod elektrodą na dwie 
części (Ryc. 1): 

— część, gdzie transport masy zachodzi głównie wskutek konwekcji (y>5); 
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Ryc. 1. Profil stężenia substancji elektroaktywnej w sąsiedztwie wirującego dysku 


— część, gdzie dyfuzja pełni główną rolę w przenoszeniu masy (y<5). 

Gdy y=5, to y=A,=1,288. Z równania (14) wynika, że graniczny prąd dysku jest 
liniową funkcją pierwiastka kwadratowego z szybkości wirowania elektrody. 
Grigory iRiddiford [18] podali bardziej dokładną postać równania (14), 
biorąc pod uwagę przy obliczeniach drugi człon w równaniu (6): 


ip, = 0,62:nF AD” : v16 2 -[1+0,3539(D/v)"* }-c°. (16) 


Ta poprawka wynosi zazwyczaj 3-5%: Ponadto wykazano [17], że różnice 
między równaniami (14) i (16) są kompensowane przez wiele czynników, m.in. 
efekty brzegowe, wpływ naturalnej konwekcji na transport masy, turbulencję, mi- 
grację jonów itp. 


3. PODSTAWY TEORETYCZNE METODY WIRUJĄCEJ ELEKTRODY 
DYSK-PIERŚCIEŃ 


Frumkin iNekrasov [19] rozwinęli koncepcję wirującej elektrody dys- 
kowej, umieszczając wokół dysku w tej samej płaszczyźnie drugą elektrodę o 
kształcie pierścienia, oddzieloną od dysku warstwą izolatora (Ryc. 2). Elektroda 
dysk-pierścień składa się z elektrody wewnętrznej (dysku) o promieniu 7,; obszaru 
warstwy izolacyjnej oddzielającej dysk od pierścienia ograniczony promieniami 


124 Jarosław Penar 


Ryc. 2. Schemat elektrody dysk-pierścień 


r, ir, elektrody zewnętrznej (pierścienia) o promieniu wewnętrznym r, i zewnę- 
trznym r; warstwy materiału izolacyjnego. RRDE charakteryzują dwa parametry 
teoretyczne: efektywność zbierania (N,) oraz efektywność ekranowania (S). Jeśli 
reakcja 


A + n,e >B (17) 
zachodzi na powierzchni dysku (przy czym początkowe stężenie substancji B jest 
równe zero) przy stałym prądzie i, , to produkt reakcji B jest przenoszony drogą 


konwekcji w głąb roztworu, a pewna jego część dociera do pierścienia. Jeśli po- 
tencjał pierścienia jest taki, że cała ilość B docierająca do niego ulega reakcji 


B:n;e—A, (18) 


to prąd pierścienia ią osiąga wartość graniczną (graniczny prąd pierścienia ) i jest 
związany z prądem dysku zależnością 
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ię, = (mz ln, )* i, No (19) 


gdzie N, jest efektywnością zbierania. Współczynnik N, wskazuje, jaka część 
wytworzonego na dysku produktu reakcji (17) przedostaje się drogą transportu 
dyfuzyjno-konwekcyjnego do powierzchni pierścienia i ulega tam reakcji (18). 
Dotyczy to oczywiście trwałych, rozpuszczalnych i nie adsorbujących się produ- 
któw reakcji. Odchylenie iz od wartości granicznej iz,, odpowiadającej Ny , wska- 
zuje na udział substancji B w następczej reakcji elektrochemicznej na dysku lub w 
roztworze. Odchylenia te mogą być też spowodowane adsorpcją B na powierzchni 
dysku. 

Zależność efektywności zbierania od parametrów geometrycznych elektrody 
podali Levich ilvanov [20]. Ci autorzy przyjęli w swych rozważaniach, że 
pierścień i przerwa izolacyjna są dużo węższe niż promień dysku, tzn. r, —7, << r, 
oraz r,— r,<< r, . Dokładniejsze wyprowadzenie zależności na efektywność zbie- 
rania podali Albery iBruckenstein [21, 22] opierając się na rozwiązaniu 
równania (3) (gdzie c jest stężeniem substancji B) przy pomocy transformacji La- 
place’a oraz funkcji Airy. To równanie zostało rozwiązane przy następujących za- 
łożeniach: 

— w układzie osiągnięty jest stan stacjonarny, tzn. dc/dt=0; 

— stężenie substancji B jest niezależne od parametru Q ze względu na symetrię 

_osiową układu; 

— współczynniki dyfuzji nie zależą od stężenia; 

— dyfuzja radialna jest zaniedbywana, ze względu na jej małą wartość w po- 
równaniu z konwekcją radialną [10]. 

Równanie (3) można zatem napisać w następującej postaci: 


v, + dc/ar + v, * ðclðy = D: (3? c/ay*), (20) 
gdzie c jest stężeniem substancji B, v, iv, dane są równaniami (4) i (6) (przy 
czym brane są pod uwagę tylko pierwsze człony tych równań). Równanie (20) 
zostało rozwiązane dla następujących warunków brzegowych: 

— c=0 dla y>% 
—(8c/dr) m = 0 dla r < r,; 
(d¢/dy),.0 = -ip /nFD 


— (Ac/dy),.0 =0 dla r, <r < ry; 
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— cz(0 dla y=Qira<r< ry 


tzn. potencjał pierścienia jest taki, że prąd pierścienia osiąga wartość graniczną 
ię. 

Rozwiązanie równania (20) dla powyższych warunków umożliwia zapisanie 
następującego wyrażenia dotyczącego efektywności zbierania elektrody N, [22]: 


N, = 1—F (a/B) + B°? - [1 —F(a)] — (1+ a+ B)** {1 -F [(a/B) : (1+ a +B)]), (21) 


gdzie 
a z(r, fr)? —.1; (22) 


B=(r;/r,)* - (rlr, )*. (23) 


Wartości funkcji F zostały podane przez autorów w postaci liczbowej dla sze- 
rokiego przedziału wartości a. i B [22]. Z równania (21) wynika, że efektywność 
zbierania elektrody N, jest funkcją jedynie parametrów geometrycznych elektro- 
dy. 

W praktyce często zdarza się, że stężenie substancji biorącej udział w reakcji na 
pierścieniu jest w głębi roztworu większe od zera. Jeżeli założy się, że reakcja 
(17) zachodzi zarówno na dysku, jak na pierścieniu przy potencjale prądu grani- 
cznego, to przy otwartym obwodzie dysku (tzn. gdy i, = 0) prąd płynący przez 
pierścień będzie miał wartość ię,. Przy zamkniętym obwodzie dysku prąd pier- 
Scienia zmniejszy się, ponieważ ilość substancji A docierającej do pierścienia 
zmniejszy się o ilość przereagowaną na dysku. Tak więc iq, < ią, , gdzie ig, jest 
granicznym prądem pierścienia, gdy na dysku płynie prąd i,, aig, jest granicz- 
nym prądem pierścienia dla i, = 0. Prąd graniczny pierścienia obniża się o war- 
tość N, *i„, a zatem 


ią = ię —No* ips - (24) 


Wzór wyrażający efektywność ekranowania (definiowaną jako wzajemny sto- 
sunek powyższych prądów pierścienia) ma następującą postać: 


S= ig, lig’, (25) 


Z rozważań hydrodynamicznych wynika, że: 
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ig, = B® ip, - (26) 


Wstawiając równania (24) i (26) do równania (25) otrzymamy: 


S=1-N,:B%. (27) 

Ponieważ N,B?*< 1, parametr S jest zawsze większy od zera. Z wyrażenia (27) 

wynika, że efektywność ekranowania przez dysk zależy wyłącznie od parametrów 

geometrycznych elektrody. Najmniejszą wartość S będą miały elektrody z wąskim 

pierścieniem i wąską przerwą izolacyjną, gdyż wtedy N, +8. Ma to duże zna- 

czenie analityczne, gdy zachodzi konieczność eliminowania wpływu substancji 
obecnej w roztworze na zachowanie się elektrody. 


4. WYKORZYSTANIE RDE DO WYZNACZANIA PARAMETRÓW 
KINETYCZNYCH I WYJAŚNIANIA MECHANIZMU REAKCJI 
ELEKTRODOWYCH 


4.1. WYZNACZANIE STAŁYCH SZYBKOŚCI REAKCJI ELEKTRODOWYCH 


W przypadku, gdy na dysku przebiega reakcja pierwszorzędowa 


Ox + ne > Red, 
a przebieg tej reakcji nie jest zbyt szybki (proces elektrochemiczny jest kontrolo- 
wany zarówno szybkością transportu, jak szybkością przeniesienia ładunku), prąd 
dysku i, (katodowy) jest określony równaniem Levicha (11): 

i, =F A*K COMA, ker C*:D—), (11a) 
gdzie k,, jest stałą szybkością katodowej reakcji elektrodowej. Równanie (11a) 
można przedstawić w następującej formie: 

i,=n*F'A' CZ: D./(Dik, + 1,61 *) Do,» w+ v") = 


=n-F-A*CQ *D, Dok, + Son) (28) 


128 Jarosław Penar 
lub 
Viy= 1/(n- F +A - kp * C8) + 1,61-v" (n F -A Do,” - CX 0). (28a) 


Analogiczne równania można zapisać dla prądów anodowych. Równania (28) i 
(28a) można stosować jedynie dla potencjałów przekraczających potencjał rów- 
nowagi E”na tyle, aby wpływ reakcji przebiegającej w drugą stronę (reakcji ano- 
dowej) był pomijalnie mały, tzn. n >>RT/anF (gdzie n=E—E ”). Przy wyznaczaniu 
k, można stosować zarówno równanie (28), jak (28a). Zastosowanie równania 
(28) zostało szczegółowo opisane przez Galusa iAdamsa [23]. Należy wy- 
znaczyć zależności prądu dysku od w'? dla różnych potencjałów. Gdy potencjał 
elektrody jest stosunkowo bliski potencjałowi formalnemu E$, natężenie prądu 
dysku jest praktycznie niezależne od w. Oznacza to, że proces kontrolowany jest 
szybkością przeniesienia ładunku, czyli wartość k, jest mała (a więc Do, /kp 
>> 0,, ). Równanie (28) można zapisać w następujący sposób: 


ip =n F- A kp CY. (29) 


Jeśli zastosujemy zależność k,, = k, ` exp[a: n -F : (E —E; )/RT], gdzie k, jest 
standardową stałą szybkością reakcji elektrodowej, a a jest współczynnikiem 
przejścia, otrzymamy 


log i, = log(n -F -A * Co, 'k,)- (a 'n'F/(2,303 -R - T) ) (E-E%). (30) 


Zależność powyższa jest prostoliniowa. Ekstrapolacja linii prostej do E= 
E? umożliwia wyznaczenie wartości k,. Z wartości kąta nachylenia prostej otrzy- 
mujemy wartość parametru œ. W podobny sposób można przeanalizować prądy 
anodowe. Aby w takiej analizie można było wykorzystywać także prądy, które za- 
leżą od w, należy stałe k, ik,, obliczać z równania (28). Wykreślając zależność 
logarytmów wyznaczonych stałych od potencjału, uzyskuje się dwie proste, przy 
czym rzędna punktu przecięcia się tych prostych przy potencjale E? jest wartością 
logarytmu stałej k, 

Wykorzystanie równania (28a) polega na wykreśleniu zależności 1/i, vs. 1/01? 
dla danego potencjału. Otrzymane proste ekstrapolowane do 1/01? 0 (tzn. w + 
0) wyznaczają wartości k, przy danym potencjale. Wartości k, i a otrzymamy, jak 
poprzednio, z analizy zależności logk, i logk,, od E. Przekształcając równanie 
(28a) do postaci: 
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1 1 1 
1/1, = GFA EC $ | aeg } (31) 


można otrzymać wyrażenie bezpośrednio określające wartość k„, mianowicie: 


ky = in, (n+ F -A - Co.) * ifipa ip) (32) 


Analogiczne równanie można zapisać dla stałej k,,. Wyrażenie podobne do 
równania (31) zostało wyprowadzone w odmienny sposób przez Małyszkę 
[24, 25]. 

Jeśli rozpatrywana reakcja elektrodowa przebiega w obydwu kierunkach, a w 
roztworze obecne są obie formy układu redoks (czyli forma Ox obok formy Red), 
to 


i,=n*F*A:D,(CQ,-C,(0))/6,,= n-F:A*D;; (Cog (0) —Cipes)/ Oper 683) 


Jeśli i , oraz ig, oznaczają katodowy i anodowy graniczny prąd dysku, to 


ię, =n- FA: DY Coo. (34a) 
i =n- F -A -Dpi ChafÒgea (34b) 


Z równań (33) oraz (34a, b) wynika, że 


Cc(0)/ C2, = (ing — ip th, = 1- iiy 


(35) 
Cz.(0)/ Choa = (i, + i) ih = 1+ Łyjiż, - 
Wyrażenie opisujące wypadkowy prąd dysku ma następującą postać: 
i,=n*F:A:(k,' Co(0)-k,' Cpe), (36) 


gdzie k,„=k,'exp[(1-a)-n'F-(E- E,)/RT] jest stałą szybkości anodowej reakcji ele- 
ktrodowej. Wstawiając zależność (35) do równania (36) otrzymamy 
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ip = n+ F +A - [hy Cy, *(l-iyfig,)-ky: Cha (1+ infig) (37) 
Ponieważ [26] 
Ki” Cor = kw" Crear (38) 


gdzie k, ik, są stałymi szybkości reakcji elektrodowej przy potencjale E”, to 


E'= E?+(R:T/(n*F))* ln (ką, kę). (39) 

Z zależności opisującej k, i ką wynika, że E = E? + (RT/nF): ln (k„/k,„). Stąd 
mamy 

n=E-E'=(RT/nF) : Intk,, ka, / (ki: Ki} (40) 


Wstawiając równania (40), (38) oraz (35) do równania (37) otrzymamy 


ip = NFA: ky Cos’ [ibs — bib, — (ih + tp) ip, * exp (nFn/RT)] . (41) 


To równanie można przedstawić w postaci 


-og(nFAk„CQ,)= log[(1/i, — 1/i5 ,) — (1/i, + Vip, ): exp (nFAn/RT)]. (41a) 


Podobne wyrażenia otrzymuje się dla stałej k,. Jeśli nFn/RT>>1 lub 
nFn/RT<<1, to otrzymane zależności przyjmują postać równania Levicha (32). 
Powyższa metoda może być stosowana jedynie wtedy, gdy rejestrowane krzywe 
i-E wykazują odpowiednie odchylenie od krzywych odwracalnych. 

Innym sposobem wyznaczania parametrów reakcji elektrodowej jest metoda 
opracowana przez Jahna i Vielsticha [27]. Stosując wzory (36) i (34a) i 
(34b) do obliczenia stężenia form Ox i Red na elektrodzie, otrzymamy następują- 
cą zależność: 


C,(0) = C2, = ip Ag, Nn :F*A: Do) > (42a) 
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oraz 


Cz(0) = Cha + i * Az (READ pn) ip, (42b) 


gdzie Ao, e= 1, 62 * Dra’ V" Wprowadzając wyliczenia dotyczące powy- 
ższych stężeń do równania (36), otrzymamy 


iz n+ F- A (kp C8, —k,, © CDNA (ky AJD): OP (Ky “Age/Dp..) 0] (43a) 


lub 


Vi, = [nFA(k, Coy — Kin Creal H+(k, Ao/D on + Kay’ Area Drea) 0712]. (43b) 


Wykreślając zależność 1/i, od w? dla stałej wartości nadnapięcia można wy- 
znaczyć wartości k,, i k,, z kąta nachylenia powyższej prostej oraz z jej ekstrapo- 
lacji do w— œ (Ryc. 3). Mając wyznaczone wartości stałych szybkości reakcji 
anodowej i katodowej dla kilku potencjałów można obliczyć standardową stałą 
szybkości reakcji elektrodowej k, (z zależności Ink, — lub Ink,, — od E), a także 


m EE Em wj | 
vfo 


Ryc. 3. Schemat zależności 1/ip vs 1/01? zgodny z równaniem (43b) , określony dla różnych 
wartości stężenia i nadnapięcia; c; > c, 1; > a 
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wartość prądu wymiany Į oraz wartość parametru a z zależności 


nFAk, = ©: exp(-anFn/RT), (44a) 


nFAk,, = i exp[(1-0)nFn/RT] . (44b) 

Powyższa metoda może być łatwo zaadaptowana do sytuacji, gdy w roztworze 
znajduje się tylko jedna z form układu redoks, np, forma Ox. Mamy wtedy: 

1fi, = [FAk,, Col" [1+ (kpAo/D oo] (45) 

To równanie po odpowiednim przekształceniu przechodzi w równanie Levicha: 


i = 1/n*F-A*k, * C2) + Vip, (31) 


Gdy w układzie badanym znajduje się tylko jedna z form układu redoks wygod- 
nie jest korzystać z metody opartej na analizie quasi-odwracalnych krzywych 
woltamperometrycznych, opracowanej przez Kutę iYeagera [28]. Ta meto- 
da opiera się na, podanej przez Korytę [29, 30], analizie równań quasi-odwra- 
calnej fali polarograficznej. Dla gałęzi katodowej to równanie ma następującą po- 
stać: 


(ip. — tp Vip * Kred Sor = €xp[nF(E —E} )/RT] + (Kz/nFk,) - exp[anF(E—E/) , (46) 
gdzie k = 0,62nFAD?*v -4/5w!2, Oznaczamy przez i, prąd odpowiadający potencja- 


łowi półfali E;, hipotetycznej odwracalnej krzywej redukcji. Jeśli posłużymy się 
zależnością 


E;, = E} + (RTF) « In(z./Ko), (47) 


to wówczas otrzymamy równanie 


k, = 0.62 © V4 - wi? - (Dg)? - (Da)? [i’ ioa). (48) 
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Sji S współczynniki kierunkowe asymptot | i 2 


Ryc. 4. Schemat analizy logarytmicznej krzywej polaryzacyjnej procesu quasi-odwracalnego na 
RDE. Krzywa 3 — krzywa doświadczalna; proste 1 i 2 — asymptoty do krzywej doświadczalnej 


Sposób wyznaczania wartości i”, oraz a wynika z ryciny 4, gdzie krzywe pola- 
rograficzne są analizowane w układzie log[ży/(i, ,— i) |vs. E. Zależność powyższa 
nie jest prostoliniowa dla procesów quasi-odwracalnych. Asymptota 1 o nachyle- 
niu charakterystycznym dla procesu odwracalnego zbiega sięz analizowaną krzy- 
wą przy potencjałach mniej ujemnych. Ta asymptota dla logiy/(i,,—i)= 0 wyzna- 
cza wartości E’,, oraz Į. Asymptota 2 jest opisana zależnością 


E = Æ + (RT/anF)- ln [ip Kip, —i,)] (49) 


i zbiega się z analizowaną krzywą przy potencjałach skrajnie ujemnych. Jej na- 
chylenie jest zależne od wartości an. Mając wyznaczone parametry i'p oraz a 
otrzymujemy wartość k, z równania (48). Ta metoda jest przydatna, gdy analizo- 
wana reakcja elektrodowa nie przebiega zbyt szybko. Jeżeli reakcja elektrodowa 
zachodzi zbyt szybko, określenie stopnia nieodwracalności nie jest precyzyjne, a 
wartość stałej k, obarczona jest dużym błędem. Z drugiej strony, gdy układ wyka- 
zuje zbyt duży stopień nieodwracalności, właściwe przeprowadzenie asymptoty 1 
jest dość trudne, co prowadzi do niezbyt precyzyjnego wyznaczenia wartości k,. 
Opisana metoda jest prosta, jednakże wydaje się, że bardziej dokładne wyniki 
uzyskuje się korzystając z zależności logk,, (lub k,,,) od E (lub od n); przeprowa- 
dzenie asymptot do krzywych eksperymentalnych obarczone jest zazwyczaj zbyt 
dużym błędem. 
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Inną metodę wyznaczania parametrów reakcji elektrodowej, opartą na analizie 
jednej (anodowej lub katodowej) gałęzi fali polarograficznej, podali Jordan i 
Javick [31]. Stosując zależność i,=n-F-A'k,' {co,(0) - exp[-a-n-F-(E- 
E})/RT — cz.„;(0) * exp[(1-0) : n: F- (E —E; /RT]) oraz równanie (50) 


ip =n: FA ` Mos* (Coe Co,(0)) = n -F -A mę: (Cę.(0)— Crea) » 50) 


gdzie m= 0,62 - D? : v16 . wi, a także biorąc pod uwagę warunek początkowy 
Chea = O (w układzie obecna jest tylko forma utleniona Ox), możemy zapisać na- 
stępującą zależność: 


ip ,/ip = (mojk,) * exp[anF(E —E,) [RT] + (m,/m,,,) ` exp [nF (E -E,)RT]+1 (51) 


Czynnik bezwymiarowy m,,/k, określa stopień nieodwracalności układu ele- 
ktrochemicznego [31]. Gdy m,/k, + 0, to fala polarograficzna jest falą chara- 
kterystyczną dla procesów odwracalnych. Gdy m,/k, > ©, to fala przybiera po- 
stać charakterystyczną dla układów nieodwracalnych. Aby obliczyć stałą k, oraz 
a., należy równanie (51) napisać w postaci 


Mo, ky = (i, — ii, — (m/m) * exp[nF(E —9)] (52) 


lub 


0,62- D$-v5- wiZ/k, = (ip, ~ ipMip — (DoD pei)” exppnF(E- E9) (52a) 


Analizując krzywe i-E dla danej wartości œ, otrzymuje się wartości ką dla sze- 
regu potencjałów. Wartość stałej k, oraz wartość parametru a. otrzymujemy, tak 
jak we wcześniej opisanych metodach, z zależności logk,, od E. 

Małyszko zaproponował zmodyfikowaną i jednocześnie uproszczoną wer- 
sję równania (52a) [24, 325]. Aby obliczyć całkowicie odwracalną krzywą kato- 
dową, należy posłużyć się następującym równaniem: 


E = E;,+(RTInF) - ln [(i,, -i’,)/i’p), (53) 


gdzie 
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Ey, = E + (RT/nF) + ln (D,„/D)?. (54) 


Prąd i”, jest hipotetycznym odwracalnym prądem dysku, płynącym przy poten- 
cjale E,„. Biorąc pod uwagę zależności (53) i (54), można równanie (52a) prze- 
kształcić w następujący sposób: 


0,62 Do, VS- wk, = (ip, - Whip — (in, - Ii", . (55) 


Podstawiając do powyższego równania wyrażenie (14) otrzymamy 


p. (n:F-A: CG: k)=(iu"i)/h-G>EpD/i'h (55a) 


Wartość drugiego wyrazu z prawej strony równania (53a) wyznacza się przez 
ekstrapolację rosnącej części krzywej woltamperometrycznej, tzn. z analizy za- 
leżności log (in, — ip)i, vs. E. Jeżeli układ jest quasi-odwracalny, to ta zależność 
jest nieliniowa. W zakresie potencjałów bardziej dodatnich (tzn. E > +0) zależ- 
ność powyższa powinna przejść w prostą opisaną równaniem (53). W tych warun- 
kach i, >i”p, a zatem asymptota do krzywej quasi-odwracalnej w tym zakresie po- 
tencjałów powinna odpowiadać procesowi całkowicie dyfuzyjnemu. Asymptota 
w punkcie przecięcia się z osią potencjału wyznacza potencjał E,„ (Ryc. 5). Z wy- 
kresu na rycinie 5 można więc odczytać wartości wyrażeń (i, ,— —i)i' (punkty 


login’ (ipr ip) 


Ryc. 5. Analiza logarytmiczna gałęzi katodowej krzywej polaryzacyjnej na RDE. Krzywa b — 
krzywa doświadczalna; krzywa a — krzywa teoretyczna zgodna z równaniem (49) 
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na asymptocie a) oraz (ip, — ż,)/i, (z krzywej b) dla danego potencjału. Znając te 
wartości można obliczyć wartości stałych ką dla różnych potencjałów, a następnie 
wartości k, oraz œ. W przypadku znacznej nieodwracalności procesu można zasto- 
sować zamiast równania (55a) równanie 


ka= lipa] Mn F A+ Co) ifn, - tp)» (32) 


ponieważ dla dużych nadnapięć 7’, = ip}. Równania analogiczne do równań (32) i 
(55a) można wprowadzić także dla procesu anodowego. Równanie ip= i}, ' {co,(0) 
/ chy exp[—a * n * F - n/RT|]-cz.(0)/cpea * xp[(1-0) * n*F - n/RT]) można prze- 
kształcić w ten sposób, by otrzymać równanie odpowiednie do wyznaczania para- 
metrów kinetycznych reakcji elektrodowej, podstawiając w miejsce stosunku 
c(0)/c® równanie (35). Otrzymamy wówczas 


1/i, = b*/{(1—-b)i9,] + (1/DŹ « 0 + IDR} + 1,61 -v - w*2/[(1 — b)nF]. (56) 
Posłużywszy się tą zależnością możemy wyznaczyć prąd wymiany oraz współ- 
czynnik przejścia (wykreślając zależność 1/i, vs. w"? dla kilku potencjałów). 
Warto na koniec omówić jeszcze stosunkowo prostą metodę zaproponowaną 
przez M iillera iPrumke [32], w której nie przyjmuje się żadnego założenia 
co do rzędowości procesu elektrodowego. Prąd dysku można ogólnie przedstawić 
następującym równaniem (analogicznym do równania (29)): 
ip=n-F-A-k, *|[ca(0)]*, (29a) 
gdzie k oznacza rzędowość reakcji elektrodowej. Wstawiając równanie (29a) do 
równania (50), będącego ogólnym wyrażeniem na natężenie prądu dysku procesu 
kontrolowanego dyfuzją, otrzymamy 
ip=n-F -A+ky* (Co): [1- i/nFAm, CO). (29b) 


Powyższe równanie można przedstawić w prostszej postaci 


ip =n*F-A*k, (eo) [1 —iplip, (29c) 
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lub 


log(i,) = login - F -A ` kp’ (Co)*] + k- log(1 —i,/i,,), (29d) 


gdzie pierwszy człon po prawej stronie równania wyraża wartość prądu aktywa- 
cyjnego. Prąd aktywacyjny, a co za tym idzie stałą k,, wyznacza się na podstawie 
przecięcia liniowych zależności log(ip)= f(1 —ip/ip,), a z nachylenia tych prostych 
można obliczyć rząd reakcji elektrodowej. 

Istnieje zatem kilka metod wyznaczania parametrów kinetycznych prostych 
procesów elektrodowych. W zależności od konkretnego przypadku można dobrać 
najbardziej wygodną i optymalną metodę. Na przykład jeśli w układzie znajdują 
się obie formy układu redoks o kontrolowanych stężeniach, a proces elektrodowy 
zachodzi stosunkowo szybko, należy posłużyć się metodą Randlesa [26] lub 
Jahna-Vielsticha [27]. Gdy układ wykazuje znaczny stopień nieodwracal- 
ności, wygodne jest stosowanie metody opartej na równaniu Levicha lub uprosz- 
czonej wersji metody M ały szki [24]. Gdy, z kolei, w roztworze występuje tyl- 
ko jedna z form układu redoks lub jedna z form tego układu jest tak nietrwała, że 
nie można kontrolować jej stężenia, należy wykorzystać metodę J ordana—Ja- 
vicka [31], Kuty-Yeagera [28]lub Matyszki [24]. 


4.2. WYZNACZANIE WARTOŚCI WSPÓŁCZYNNIKA PRZEJŚCIA 
Ogólny współczynnik przejścia wyznacza się na podstawie analizy krzywej i-E 


w układzie współrzędnych log[i,/i,,— ż,)] vs. E. Ponieważ równanie nieodwracal- 
nej krzywej i,—E ma postać 


E + E? + (RT/anF) - In (D/k,8) + (RT/anF) An [i,/(i,, -i,)], (57) 


gdzie 


E? + (RT/anF) - In(D/k,d) = Ey» (58) 


zatem współczynnik nachylenia prostej przedstawiającej zależność prąd-poten- 
cjał w powyższym układzie współrzędnych powinien być równy 0,059/anF V na 
jednostkę logarytmu. Także liniowa zależność potencjału półfali E,, od log(w) po- 
zwala na wyznaczenie wartości a. Współczynnik przejścia można także obliczyć 
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postuzywszy się zależnością logky, (lub logk,,,) od potencjału w sposób omówiony 
już w poprzednim rozdziale. 
4.3 WYZNACZANIE POTENCJAŁU FORMALNEGO 

Potencjał formalny nie jest co prawda wielkością kinetyczną, jednakże znajo- 
mość jego wartości jest niezbędna przy wyznaczaniu standardowej stałej szybko- 
ści reakcji. Potencjał formalny można wyznaczyć przy pomocy RDE na podsta- 
wie analizy natężenia prądu dysku w zależności od potencjału [33] w układzie log 
iy/(in— ip) vs. E. Należy najpierw wyznaczyć wartość odwracalnego potencjału 
półfali E;,, a następnie stosować równanie (54). Gdy nachylenie prostej wyzna- 
czonej przy pomocy powyższej zależności jest różne od wartości teoretycznej 
2,303RT/nF, może to świadczyć o nieodwracalności badanego procesu i potencjał 
formalny należy wyznaczyć inną metodą. 

4.4. OKREŚLENIE STOPNIA NIEODWRACALNOŚCI PROCESU 


Równanie Levicha (31) można zapisać w postaci 


infin, = 5k, /(5 * kyt Do,)- (59) 
Proces będzie odwracalny, gdy spełniona zostanie zależność 


Ò- ky>>Do, . (60) 

Zakładając wartość Dp, równą 105cm*/s, mamy ô: k,, > 10*%cm*/s. Proces ele- 
ktrodowy będzie kontrolowany wyłącznie szybkością przeniesienia ładunku, jeśli 
D.>>ky * 5. (61) 


W przypadku procesu nieodwracalnego mamy 


infin, = ky * SID. (62) 


Gdy D,,= 10°cm?/s, to 5k, < 10cm?/s. Na rycinie 6 przedstawiono zależność 
natężenia prądu dysku od œ? dla procesu, który dla małych wartości w jest od- 
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Rys. 6. Zależność prądu dysku od w!?. Potencjał elektrody leży w obszarze formowania się fali. Pro- 

ces elektrodowy jest quasi-odwracalny. Strefa I — reakcja elektrodowa kontrolowana szybkością 

transportu; strefa II — reakcja elektrodowa kontrolowana zarówno szybkością transportu, jak i szyb- 

kością przeniesienia ładunku; strefa III — reakcja elektrodowa kontrolowana szybkością 
przeniesienia ładunku 


wracalny (strefa I), przy szybkościach średnich staje się quasi-odwracalny (strefa 

II), a przy szybkościach bardzo dużych staje się nieodwracalny (strefa III). 
Trochę inaczej warunek nieodwracalności został sformułowany przez J o r d a - 

na [34]. Równanie (51) można, przy założeniu że mo, = Mg» zapisać w postaci 


i, = i, (mak, + 2). (63) 


Prąd jest kontrolowany szybkością transportu masy, gdy m,,/k,<<2 (tzn. gdy 
D>>2k,o,). Gdy m,/k,;>>2 (tzn. D<< 2k,0,,), to prąd kontrolowany jest szybko- 
ścią przeniesienia ładunku, a więc proces jest nieodwracalny. Nieodwracalność 
procesu elektrodowego powoduje, że potencjał półfali zmienia sięwraz ze zmianą 
szybkości wirowania elektrody. Prócz tego zależność logi,/(i,, — ip) vs. E ma na- 
chylenie różne od wartości teoretycznej. Określenie stopnia nieodwracalności 
procesu elektrodowego ma duże znaczenie, jeśli badamy kinetykę reakcji elektro- 
dowej. Pomiary stałych szybkości reakcji można przeprowadzić tylko w warun- 
kach procesu quasi-odwracalnego. Należy zatem tak dobrać warunki prowadzenia 
pomiarów, aby badany proces miał charakter quasi-odwracalny. Z kolei wyzna- 
czanie wartości potencjału formalnego można przeprowadzić tylko wtedy, gdy 
proces jest odwracalny. Zmieniając odpowiednio szybkość wirowania elektrody 
można w dość dużym zakresie zmieniać parametr 6, a przez to wpływać na chara- 
kter badanego procesu. 
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4.5. WYKORZYSTANIE RDE DO BADANIA REAKCJI HOMOGENICZNYCH 


Ważnym zastosowaniem RDE jest badanie kinetyki reakcji chemicznych za- 
chodzących w głębi roztworu. Właściwości układu badawczego z RDE pozwalają 
na rozwiązanie równań transportu nawet w tym przypadku, gdy w układzie wy- 
stępują skomplikowane reakcje chemiczne. 

Wówczas, gdy substancja biorąca udział w reakcji elektrodowej wchodzi także 
w równowagi chemiczne w głębi roztworu, równanie dyfuzyjno-konwekcyjne (8) 
należy uzupełnić o kinetyczny człon opisujący powyższe równowagi. Jeśli rozpa- 
trzy się pojedynczą reakcję chemiczną, w której biorą udział substancje A,, A,, ... 
A,, to jej równanie chemiczne można zapisać w następujący sposób: 


a v, A, =0, (64) 


gdzie v, oznacza odpowiedni współczynnik stechiometryczny. Stała równowagi 
reakcji wyraża się równaniem 


o= H ov,- (65) 


Gdy któraś z substancji A,, .... A, bierze udział w reakcji elektrodowej, równo- 
waga reakcji chemicznej — dana równaniem (65) — w pobliżu elektrody zostaje za- 
chwiana, a odchylenie od równowagi można wyrazić funkcją [11] 


Pf (Cy C) =P „dh c(0)v, — r: a 14,0) Vp (66) 


gdzie p jest stałą szybkości reakcji chemicznej, a c,(0) oznacza stężenie substancji 
A, w pobliżu elektrody. 

Jeśli reakcji elektrodowej ulega substancja A,, to reakcja chemiczna w głębi 
roztworu powoduje wzrost lub zmniejszenie stężenia tej substancji w pobliżu ele- 
ktrody, a co za tym idzie — zmieniają się warunki transportu dyfuzyjno-konwe- 
kcyjnego. Dla procesu stacjonarnego równanie dotyczące transportu można zatem 
zapisać w następującej postaci: 


D, @efdy*) = v, * (de,/dy) + Q,, (67) 
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gdzie 


Q= 'P'/(--.C,). (68) 


Warunki brzegowe są następujące: c, = c? dla y + © oraz c, = 0 dla y= 0, gdy 
proces elektrodowy przebiega odpowiednio szybko. Ponieważ substancje A,, , nie 
biorą udziału w procesie elektrodowym, to 


(de/dy)„„= (dejdy),„= ... = (de,/dy) „, (69) 

Metody rozwiązywania równania (67) dla różnych typów reakcji chemicznych 
biorących udział w procesie elektrodowym zostały opisane w monografii (11). 
Poniżej przedstawiono kilka przykładów częściej spotykanych w praktyce do- 
świadczalnej przypadków, które ilustrują przydatność RDE w tego typu bada- 
niach. 


4.5.1. BADANIE KINETYKI REAKCJI CHEMICZNEJ POPRZEDZAJĄCEJ REAKCJĘ ELEKTRODOWĄ (PROCES TYPU CE) 


Proces typu EC można przedstawić równaniem 


A, $>A4, 


przy czym (70) 
A,+ ne — produkty. 


Układ równań opisujących stacjonarny transport masy i odpowiednie warunki 
brzegowe jest następujący (przy założeniu D, = D,): 


D : (cjd?) = v, (dc,/dy) — p * (c, — pc,), (71) 
D: (#c,Jdy*) = v,* (dc./dy) + p: (c, — pc,), (72) 
c,=Q;c, =e dla y > ©, (73) 

c, =0; (dc./dy), „= 0. (74) 


Rozwiązanie równań (71-72) zostało podane przez Koutecky ego iLe- 
vicha (41) w takiej postaci: 
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Ryc. 7. Schemat zależności ip ,/w'” vs ip, zgodnie z równaniem (79) dla 6, = 0 i 6,> 0 


j= D+ (+ MoM - [DKU +0)))2+ 1] *6)= 
| (75) 
= jool: (DK + o) 1), 


gdzie jp oznacza przepływ masy w kierunku elektrody, a jp, jest graniczną warto- 
ścią jp, wielkość (D/p(1 + o)]*? ma wymiar [cm] i traktowana jest jako grubość 
warstwy kinetycznej 5,. Rozwiązanie powyższe zostało wprowadzone przy zało- 
zeniu, że reakcja chemiczna przebiega odpowiednio szybko, tzn. 0 < ô, < 5. Rów- 
nanie (75) można przedstawić w postaci następującej zależności między natęże- 
niem prądu kinetycznego i granicznego prądu dla substancji A, lub A,: 


ip, = (ip) * S/[1 + 08/8]. (76) 


To równanie jest słuszne dla o>>1, tzn. dla pomijalnie małych stężeń substancji 
A, w roztworze. Gdy o< 1, efekty kinetyczne są bardzo małe i prąd osiąga wartość 
prądu granicznego (i,,), [11]. Prąd kinetyczny opisany równaniami (75, 76) jest 
mniejszy niż graniczny prąd dyfuzyjny substancji A,. Zależność pomiędzy i, , oraz 
(i,) wyraża równanie 


ing = (ip,)-/(1 + O * 8/8). : (7) 
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Dla praktycznych zastosowań dobrze jest równanie (77) przedstawić w takiej 
postaci: 


ipy/ (inz) a 


ip (ina 


a:0/0= flo”). (78) 


To równanie pozwala wyznaczyć wartość 6,, jeśli znana jest stała równowaga 
chemicznej o. Powyższe równanie można także zapisać w postaci następującej: 


i, Jo? =B+ C* ipl? (79) 
gdzie B = 0,62: n: F-Ds:v6. Gl, c =p"? - D/(1,61 - v!5). 


Zależność i, ,/w!? = flip,) jest zależnością prostoliniową i wartość p można wy- 
znaczyć z nachylenia prostej. Do posłużenia się równaniami (78,79) nie musimy 
mieć danych dotyczących stężenia substancji A. Równanie (79) ma także funkcję 
diagnostyczną. Jeśli w badanym procesie elektrochemicznym występuje równo- 
waga chemiczna poprzedzająca reakcję elektrodową, przy czym reakcja przebiega 
tak szybko, że spełnione są warunki: > 0 oraz 0,< ô, to wykres zależności i/w'” 
vs. i jest linią prostą równoległą do osi rzędnych. Gdy reakcja chemiczna przebie- 
ga bardzo powoli, to >> ð, a człon kinetyczny w równaniach (71-72) jest mały, 
oznacza to, że reakcja chemiczna nieznacznie tylko wpływa na wielkość granicz- 
nego prąd dysku. 

W powyższych równaniach przyjęto, że współczynniki dyfuzji substancji bio- 
rących udział w procesie elektrochemicznym są jednakowe. Dogonadze (42) 
rozważał bardziej ogólny przypadek procesu (70), w którym współczynniki dyfu- 
zji substancji A, i A, są różne, przy czym, jak poprzednio, podane rozwiązanie do- 
tyczy reakcji chemicznej przebiegającej szybko, tzn. d,<<0. Kirjanow iFili- 
nowskii (43) podali przybliżone rozwiązanie równań (71-72) bez przyjmowa- 
nia ograniczeń co do szybkości przebiegu reakcji chemicznej. Otrzymane przez 
nich wyrażenie dotyczące kinetycznego prądu dysku ma następującą postać: 


(i,,),(1 + 0) 


= —— L (80) 


o[3,1 + $0(1 + AIA’ 
me 1,65 + 87%0(1 + db 


ok 


W przypadku reakcji chemicznych przebiegających szybko (tzn. >> D/[p(1 + 
0)]) równanie (80) przechodzi w równanie (77) z dokładnością 6%. Jako reakcja 
poprzedzająca reakcję elektrodową może występować także reakcja dimeryzacji 
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A, £ 2A,, ° 
przy czym (81) 
A, + ne — produkty. 


Układ równań dla tego procesu ma postać 


D(#c,/dy?) = v,(dey/dy) — 2p(c, — pe), 


D(Pejdy?) = v,(dcyjdy) + p(z- på). 


Warunki brzegowe tego procesu są takie same, jak poprzednio (równania 73— 
74). Jeśli reakcja chemiczna przebiega szybko, tzn. gdy ô, << ô (czyli cż>> c9), 
wówczas rozwiązaniem tego układu równań jest następująca zależność [11]: 


1- (ip g ci 


ing Ain) P “BY, © 


gdzie ô, = [D/2poc? ]"?, a (ip). = ZnFDcy6. W tym przypadku zależność pomię- 
dzy prądem a szybkością wirowania elektrody jest bardziej złożona niż w przy- 
padku reakcji jednoczasteczkowej. Do gonadze (42) rozważał ten przypadek, 
gdy współczynniki dyfuzji substancji A, iA, nie były sobie równe. 


4.5.2, BADANIE KINETYKI RBAKCJI CHEMICZNEJ PRZEBIEGAJĄCEJ RÓWNOLEGLE Z REAKCJĄ ELEKTRODOWĄ (REAKCJE 
KATALITYCZNE) 


Proces katalityczny można przedstawić w postaci 


A, £ne > A, 
A, $$, A; (84) 


Powyższy proces opisywany jest równaniami (71-72), a warunki brzegowe na 
dysku są następujące: 
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c, =0; (dc,/dy),_,= (dc,/dy),_, dla y = 0. (85) 


Przy założeniu, że reakcja chemiczna przebiega odpowiednio szybko, rozwią- 
zanie powyższego układu ma postać (39): 


ing =n “F(R +) [PDK + OJJ. (86) 


Z równania (86) wynika, że prąd katalityczny, po pierwsze — nie zależy od 
szybkości wirowania elektrody, po drugie — określony jest wielkością stałej szyb- 
kości reakcji przemiany produktu reakcji elektrodowej A, w substancj¢A,. Nieza- 
leżność prądu od w wynika z tego, że 6, << ô (tzn. reakcja chemiczna przebiega 
bardzo szybko), co oznacza, że reakcja chemiczna przebiega blisko powierzchni 
elektrody, w strefie, gdzie transport konwekcyjny jest niewielki. 

Ten układ równań można rozwiązać nie przyjmując warunków co do szybkości 
reakcji chemicznej, przy czym to rozwiązanie jest rozwiązaniem przybliżonym 
[43]. Wyrażenie dotyczące prądu katalitycznego ma w tym przypadku postać: 


Ryc. 8. Proces katalityczny. Schemat zależności prądu dysku od w". Krzywa 1 — szybka reakcja 
katalityczna (84) zgodna z równaniem (86); krzywa 2 — zgodna z równaniem (87); krzywa 3 — w 
układzie nie występuje reakcja katalityczna 
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 _moą L65+p(1+0)6b 
pk= (i)! [3,1 +P(1+0)00]*7 (87) 


Z powyższego wyrażenia wynika, że dla reakcji chemicznych przebiegających 
szybko — tj. gdy p(1 + 0)8?/D>> 1) — równanie (87) przechodzi w równanie (86). 
Przy reakcjach chemicznych przebiegających powoli — tj. gdy p(1 + 0)07/D<< 1) 
— równanie (87) przechodzi w wyrażenie dotyczące prądu granicznego dysku (i, ,), 
(z dokładnością do 6%). Oznacza to także, że 

1. dla małych szybkości wirowania elektrody, gdy graniczna warstwa dyfuzyj- 
na jest duża i spełniony jest warunek 6 >> ô, prąd katalityczny nie zależy od w; 

2. dla dużych szybkości wirowania elektrody, gdy grubość warstwy dyfuzyjnej 
znacznie maleje i zachodzi zależność ô, >> 0, prąd katalityczny zaczyna rosnąć 
proporcjonalnie do pierwiastka kwadratowego z szybkości rotacji elektrody. 

Na rycinie 8 przedstawiono schematycznie zależność i,, od œ? dla trzech 
przypadków ilustrujących powyższe konkluzje. Wykres 1 pokazuje zależność i,, 
od w'? według równania (86), wykres 2 zależność i,, od w? według równania 
(87), a wykres 3 zależność prądu dysku od w** (w przypadku, gdy reakcja chemi- 
czna w układzie nie występuje). Z początkowej części wykresu 2 i zależności (86) 
można obliczyć wartość stałej p. 

Haberland i Landsberg (45) rozwiązali równanie transportu dla roz- 
patrywanego przypadku, przy założeniu, że przebieg reakcji katalitycznej nie jest 
zbyt szybki. Otrzymane przez nich równanie opisuje zależność prądu od w w for- 
mie wygodnej w praktyce: 


ilin), = [8p"7/D"7] « cotanh (ò - p'/D") = x - cotanh (x), (88) 


gdzie x = 1,61 (v/D)"* (p/w)'”, przy czym 


x= [ip /(tp,),] * tanh (x). (89) 


Należy zatem wyznaczyć wartości prądu dla (i,,), dla odpowiednich wartości 
w, gdy proces katalityczny nie zachodzi, a następnie zmierzyć wartości prądu 
kinetycznego przy tych samych warunkach w. Znając wartości i,, oraz (i,,), dla 
różnych w, oblicza się wartość x z równania (89), a z zależności x = f(w'?) wyzna- 
cza się wartość p. 

Stosunkowo prostą metodą wyznaczania parametrów kinetycznych jest metoda 
opracowana przez Bruckensteina iPragera [46], oparta na koncepcji 
warstwy dyfuzyjnej Nernsta. Wyniki otrzymane powyższą metodą są porówny- 
walne z wynikami otrzymanymi ze wzorów analitycznych. Stosując metodę Bru- 
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ckensteina i Pragera do procesu katalitycznego, otrzymuje się następującą zależ- 
ność: 
R+ x%-p-R?-1=0, (90) 


gdzie R = 0,0 = (ip4/ipu, a X = 1,61D$* v"7/D,. To równanie można przedstawić 
w postaci takiego równania prostej: 


(1-RJR* = xp. | (91) 


Łatwo teraz można obliczyć p, znając wartości (i,,), oraz ip, wyznaczone 
w sposób opisany powyżej. Gdy mamy do czynienia z reakcją drugorzędową typu 


2A, $34, ; A, £ ne A, (92) 
musimy rozwiązać następujący układ równań: 


D(@e,/dy?) = v,(de,/dy) — p(eż — oc) 
D(dc dy?) = v,(de,Jdy)+ 2p(cż —0c,), 


przy czym dane dotyczące warunków brzegowych pochodzą z równań (73) i (85). 
Rozwiązaniem takiego układu równań — przy założeniu, że 2c? c9, (tzn. stężenie 
substancji A, w głębi roztworu jest pomijalnie małe) — jest wyrażenie 


ing 2 (ipn) 


i a (3/2)*2 (2 ¢9/c9)( 8/8,), (94) 


gdzie 5,= (D/po)'”. Gdy wartość p jest duża, natężenie kinetycznego prądu dysku 
osiąga wartość maksymalną 2(ip,),. Ulstrup [47] podał analityczne rozwiąza- 
nie powyższego układu równań przy zaniedbaniu członów konwekcyjnych oraz 
przy założeniu, że równowaga reakcji chemicznej jest przesunięta całkowicie 
w prawo (c3 = 0). Ma ono następującą postać: 


1—y = 1,225 - (6/8) y, (95) 
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gdzie y = jp,(2Dci/6), jo, jest gęstością kinetycznego prądu dysku, a 6, = 
(D/kc,(0))'2= = 6*6 [D6/pcty]!?. Ta metoda pozwala na wyznaczenie stałej szyb- 
kości reakcji homogenicznej dla p = 10° dm*/mol - s, co odpowiada stosunkowi 
infin) = 2. Holub [48] rozwiązał powyższy problem numerycznie dla 
wszystkich możliwych wartości (ip, — in.) /iņņ= J w zakresie 0 < J s 1. W arty- 
kule [48] zamieszczono tabelę zawierającą wartości parametru 4 = 
[2c%p/0,51%%w] * [v/D]** w zależności od J dla x < 50. 


4.5.3. BADANIE KINETYKI REAKCJI CHEMICZNYCH PRZEBIEGAJĄCYCH PO REAKCJI ELEKTRODOWEJ (PROCES TYPU EC) 


Proces typu EC przedstawiają następujące równania: 


A,+ne <A, 
A, EA, 

Produkt reakcji elektrodowej A, ulega reakcji chemicznej w pobliżu elektrody. 
Wpływ reakcji chemicznej na kinetykę reakcji elektrodowej będzie największy 
wtedy, gdy przebieg etapu elektrodowego jest szybki (reakcja elektrodowa jest 
odwracalna), a równowaga reakcji chemicznej jest przesunięta w kierunku two- 
rzenia substancji A,, która nie bierze bezpośrednio udziału w reakcji elektrodo- 
wej (tzn. o >> 1). Układ równań różniczkowych opisujących ten proces jest taki 
sam jak układ (71-72), przy czym stężenie substancji A, spełnia zależność 


D(ac,/dy*) = v,(de,/dy), (97) 


a warunki brzegowe są następujące: 


c =c]; 0,=C3C,=C3 dla y — ©; 
(des/dy) „, = -(dc,/dy), o; (dc./dy) = 0; 
c,(0) = c,(0)exp[nF(E - EẸ)/RT] dla y —> ©. (98) 


Rozwiązaniem powyższego układu równań jest wyrażenie [11] 
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A= io Ain) | (6/5) - exp[nF(E — E9) /RT], (98) 


ing Nin: 

gdzie 5, = (D/po)'”, a E? jest potencjałem formalnym układu redoks. Z równania 

(98) wynika, że wartość potencjału półfali E,, w obecności procesu chemicznego 

jest dana zależnością 

(E;n), = E? + (RTInF)* In[(po/D)28]. (99) 

Wynika stąd, że potencjał półfali dla procesu EC zostaje przesunięty względem 

potencjału półfali dla procesu bez reakcji chemicznej, gdy E,, = Ev. Biorąc pod 

uwagę wielkość przesunięcia potencjału półfali można wyznaczyć stałą szybkości 
reakcji p. W przypadku reakcji dwucząsteczkowej, tzn. gdy zachodzi proces 

2A, FA, (100) 


rozwiązaniem układu odpowiednich równań jest zależność [29]: 


t- ing Ain) z 


Z DL 
Potencjał półfali wyraża się następującym wzorem: 


(E,),= E? + (RT/3nF) « n[po c96/3D). (102) 


Potencjał półfali zależy tu także od stężenia depolaryzatora. Z wielkości prze- 
sunięcia potencjału półfali można wyznaczyć wartość p. W przypadku procesu ty- 
pu ECE, to znaczy 

A, ż n.e >A, FA, +ne>A, — (103) 
gdy n, = n, = n, a szybkość przemiany A, jest większa niż przemiany A, prąd dys- 
ku dany jest równaniem [48] 


ing = (i); * [2 —exp(-pot MD. (104) 
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Gdy reakcja chemiczna przebiega bardzo szybko (to jest gdy po>D/5), wypad- 
kowy prąd dysku odpowiada wymianie 2n elektronów (człon eksponencjalny 
w równaniu (104) znika). Dla po = 0 wielkość prądu dysku odpowiada wymianie 
n elektronów. Dla 0< po< <D/8 obserwowana liczba elektronów wymienianych 
w procesie zmienia się od 2n do n. Tak więc równanie (104) można przedstawić 
w postaci 


in, Kinz); = n/n = 2— exp[-pod/TID], (105) 


umożliwiającej wyznaczenie wartości po. 


5. WYKORZYSTANIE RRDE DO USTALANIA MECHANIZMU I PARAMETRÓW 
REAKCJI ELEKTRODOWYCH 


Metoda RRDE jest często stosowana do badania kinetyki złożonych reakcji 
elektrodowych, w których pierwotne produkty reakcji elektrodowej ulegają dal- 
szym przemianom. Dwa niezależne układy elektrodowe umożliwiają wyznacze- 
nie parametrów kinetycznych różnych typów reakcji elektrochemicznych. Meto- 
da RRDE umożliwia wyznaczanie krzywych polaryzacyjnych na dysku i na pier- 
Scieniu, badanie ekranowania pierścienia przez dysk, wyznaczanie zależności 
wielkości granicznego prądu pierścienia od takich parametrów, jak: prąd dysku, 
potencjał dysku, szybkość rotacji elektrody, czas, stężenie depolaryzatora. Na 
podstawie takich zależności można wyznaczyć trwałość badanych reagentów 
i odróżnić produkt przejściowy od końcowego. Można na ich podstawie ustalić, 
czy proces związany z ubytkiem produktu pośredniego ma charakter chemiczny 
czy elektrochemiczny oraz gdzie zachodzi: na powierzchni elektrody czy w głębi 
roztworu. Te dane są także pomocne przy wyznaczaniu rzędowości reakcji pro- 
duktu pośredniego. 


5.1. USTALANIE MECHANIZMU REAKCJI METODĄ RRDE 


5.1.1. ZALEŻNOŚĆ EFEKTYWNOŚCI ZBIERANIA RRDE OD POTENCJAŁU DYSKU 


Kształt zależności N od E, dla różnych przypadków przedstawiono na rycinie 
9 [50]. Zmiany potencjału dysku powodują zmiany w stężeniu produktu przej- 
ściowego i końcowego w warstwie przyelektrodowej. Jeżeli na dysku powstaje 
trwały produkt, to zmiana potencjału dysku zwiększa szybkość tworzenia się tego 
produktu, lecz efektywność zbierania elektrody pozostaje stała (wykres 1 na ryci- 
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nie 9) i równa jest N,. W przypadku, gdy produkt reakcji dyskowej ulega w roz- 
tworze reakcji chemicznej I rzędu, efektywność zbierania elektrody także nie 
zmienia się, jednak jej wartość jest mniejsza od N, (wykres 2 na rycinie 9). Gdy na 
pierścieniu rejestrowanych jest kilka produktów, a jeden z nich jest produktem 
przejściowym, który przetwarzany jest na produkt końcowy na drodze następczej 
reakcji elektrochemicznej na dysku, to szybkość tej przemiany zwiększa się wraz 
ze wzrostem potencjału bardziej niż szybkość tworzenia produktu pośredniego. 
Dlatego wartość efektywności zbierania produktu przejściowego maleje (wykres 
3 na rycinie 9), a efektywność zbierania wyznaczona dla produktu końcowego 
rośnie (wykres 4 na rycinie 9). Podobne zależności wykrywa się wówczas, gdy 
przetwarzanie produktu przejściowego w końcowy odbywa się na drodze nastę- 
pczej reakcji chemicznej o wyższej rzędowości. Z przedstawionych powyżej za- 
leżności wynikają kryteria, dzięki którym można odróżnić produkt końcowy od 
przejściowego oraz wskazać chemiczny lub elektrodowy charakter przemiany 
produktu przejściowego. 


5.1.2. ZALEŻNOŚĆ EFEKTYWNOŚCI ZBIERANIA OD STĘŻENIA DEPOLARYZATORA 


Z braku zależności między, z jednej strony, wartością efektywności zbierania 
elektrody obliczonej dla produktu pośredniego i końcowego, a z drugiej — stęże- 


Ep 
Rys. 9. Schemat zależności efektywności zbierania elektrody N od potencjału dysku Ep. Krzywa 1 — 
brak reakcji w roztworze, produkt reakcji elektrodowej jest trwały; krzywa 2 — produkt reakcji ele- 
ktrodowej ulega w warstwie dyfuzyjnej reakcji homogenicznej I rzędu; krzywa 3 — produkt genero- 
wany na dysku jest nietrwały i ulega następczej reakcji elektrodowej; krzywa 4 — efektywność zbie- 
rania dla produktu końcowego reakcji elektrodowej na dysku 
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niem depolaryzatora wynika, że przemiana produktu przejściowego w końcowy 
zachodzi w reakcji pierwszego rzędu. Spadek wartości N, wyznaczonej dla produ- 
ktu końcowego, wraz ze wzrostem stężenia depolaryzatora, potwierdza tezę, że 
rząd reakcji jest większy od I. 


5.1.3. ZALEŻNOŚĆ GRANICZNEGO PRĄDU PIERŚCIENIA OD PRĄDU DYSKU 


Jeśli produkt tworzony na dysku jest trwały, to wykresem zależności i,, od i, 
jest prosta przechodząca przez początek układu o nachyleniu równym N, (wykres 
1 na rycinie 10). W zakresie wyższych prądów dysku może nastąpić odchylenie 
od tego przebiegu i spadek wartości prądu i,,, co jest związane z obniżeniem wy- 
dajności tworzenia produktu na dysku. Prąd dysku, przy którym następuje opisa- 
ne zjawisko nazywa się krytycznym prądem dysku. 

Gdy produkt generowany na dysku ulega w roztworze reakcji homogenicznej 
I rzędu, to prąd graniczny pierścienia zmienia się proporcjonalnie do zmian prądu 
dysku, przy czym wykres zależności ię, od i, ma nachylenie mniejsze niż war- 
tość N, (wykres 2 na rycinie 10). 

Gdy produkt wytwarzany na dysku ulega reakcji homogenicznej II rzędu, to za- 
leżność ią, od i, ma postać krzywej 3 na rycinie 10. Przy małych prądach dysku 
praktycznie cały generowany produkt przejściowy ulega reakcji homogenicznej, 
zanim osiągnie brzeg pierścienia; tak więc i, = 0. Ze wzrostem wartości i, ob- 


|IRL| 


| Ipi 


Ryc. 10. Zależność granicznego prądu pierścienia ip; od prądu dysku ip; krzywa a — produkt re- 

akcji dyskowej jest trwały; krzywa b — produkt generowany na dysku ulega w warstwie dyfuzyjnej 

reakcji homogenicznej I rzędu; krzywa c — produkt generowany na dysku ulega w warstwie 
dyfuzyjnej reakcji homogenicznej II rzędu 
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szar drugorzędowej reakcji zbliża się do wewnętrznej krawędzi pierścienia, gdzie 
część produktu przejściowego ulega reakcji elektrodowej, tak więc i,,> 0, przy 
czym zależność ip, od ż, jest nieliniowa. Gdy obszar reakcji przesunie się poza 
zewnętrzną krawędź pierścienia, to do pierścienia dostarczany jest produkt przej- 
ściowy, który po drodze nie ma możliwości wejścia w reakcję homogeniczną, za- 
tem wykres zależności ip, od i, w tym obszarze prądów dysku przechodzi w linię 
prostą o nachyleniu N,. 


5.1.4. ZALEŻNOŚĆ EFEKTYWNOŚCI ZBIERANIA OD SZYBKOŚCI WIROWANIA ELEKTRONÓW 


Badanie zależności N od w (przy stałych pozostałych parametrach) umożliwia 
przede wszystkim określenie, na ile stabilny jest produkt rejestrowany na pierście- 
niu. Jeśli badany produkt jest trwały, to efektywność zbierania elektrody dla tego 
produktu jest niezależna od w (wykres 1 na rycinie 11). Podobną zależność można 
otrzymać dla produktu przejściowego, którego trwałość przekracza znacznie czas 
potrzebny na pokonanie drogi z dysku do pierścienia (w zakresie wykorzystywa- 
nych szybkości wirowania elektrody). Efektywność zbierania na pierścieniu pro- 
duktów nietrwałych rośnie wraz ze wzrostem w'2, dążąc do wartości N, (krzywa 
2 na rycinie 11). Spadek efektywności zbierania wraz ze wzrostem w'” jest zazwy- 
czaj związany z pojawianiem się stabilnego, rejestrowanego na pierścieniu, pro- 
duktu wytworzonego z niestabilnych substratów generowanych na dysku (wykres 
3 na rycinie 11). Rodzaj zależności efektywności elektrody od szybkości jej wiro- 


Ryc. 11. Zależność efektywności zbierania elektrody N od w'?, Krzywa 1 — produkt reakcji elektro- 
dowej na dysku jest trwały; krzywa 2 — produkt reakcji elektrodowej jest nietrwały; krzywa 3 — 
trwały produkt reakcji elektrodowej wytwarzany jest z mało stabilnych substratów 
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wania dostarcza także informacji o rzędowości reakcji zaniku niestabilnego pro- 
duktu generowanego na dysku oraz pokazuje, czy proces ten jest hetero- czy 
homogeniczny. Liniowa zależność N* od w'1? odcinająca na osi odciętych wartość 
1/N, wskazuje, że reakcja zaniku produktu pośredniego ma charakter procesu 
heterogenicznego pierwszego rzędu (wykres A na rycinie 12). Prostoliniowa za- 
leżność N —w'! sugeruje, że ten proces jest reakcją homogeniczną pierwszego 
rzędu (wykres B na rycinie 12). 


Ryc. 12. Zależność odwrotności efektywności zbierania od 1/0*? i od 1/w. 

A-— produkt pośredni reakcji elektrodowej ulega reakcji heterogenicznej I rzędu; 

B — produkt pośredni reakcji elektrodowej ulega reakcji homogenicznej I rzędu 
W przypadku wieloetapowych reakcji elektrochemicznych zostały opracowane 
specjalne kryteria diagnostyczne, które umożliwiają określenie, czy na elektro- 
dzie dyskowej przebiega reakcja bezpośrednio prowadząca do stabilnego produ- 
ktu końcowego, czy reakcje równoległe [51]. Można taki przypadek przedstawić 

w następującej postaci: 


reakcje na dysku A PD; (107a) 
ApoB; (167b) 
BC; (107c) 
reakcja w roztworze B >E; (107d) 
reakcja na pierscieniu B >F; (107e) 


C i D są elektronieaktywne; a także może zachodzić przypadek, że C = D. Dla 
takiego procesu wyrażenie dotyczące odwrotności efektywności zbierania ele- 
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ktrody (czyli stosunek prądu dysku do granicznego prądu pierścienia) ma postać 
[51] 


UN = infig, = (x + 1YN + [(x + 2) - RIN] * Vo”, (108) 


gdzie k’ = 1,66 - D% v'*k,, a x= i/i (i, i i, są prądami związanymi odpowiednio 
z reakcją (107a) i (107b)). Czynnik x określa, czy zachodzi reakcja równoległa 
(107a). Dla reakcji bezpośredniej x = 0 (tzn. i, = 0), natomiast dla reakcji równo- 
ległej x = 1 (tzn. i, = i). Zależność 1/N od ay” jest liniowa; punkt przecięcia z osią 
odciętych określa wartość (x + 1)/N,. W opisywanym przypadku można rozpatry- 
wać pięć następujących sytuacji: 

1. występuje tylko reakcja (107a) i wtedy i, = 0; 

2. występuje tylko reakcja (107b), tak więc k, = 0 i 1/N= 1/N,. Zależność N od 
w-2 jest linią prostą równoległą do osi rzędnych (wykres A, na rycinie 13); 


SARE 


* + 
R X 
Ja a? 


Ryc. 13. Zależność 1/N od 1/w": a — reakcja elektrodowa przebiega jednotorowo (107b) z wytwo- 

rzeniem produktu pośredniego, który nie ulega dalszym przemianom; b — reakcja elektrodowa prze- 

biega jednotorowo (107b) z wytworzeniem produktu pośredniego, który ulega dalszej przemianie 

(107c); c — produkt pośredni wytwarzany jest w reakcji równoległej (107a i 107b) i nie ulega 

dalszym przemianom; d — produkt pośredni wytwarzany jest w reakcji równoległej (107a i 107b) i 
ulega dalszym przemianom (107c) 
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3. zachodzą tylko reakcje (107b) i (1076), bez reakcji równoległej (107a), czyli 
x =0. Z równania (108) wynika 


1/N = 1/N, + 2k IN, 0!2). (109) 


Tak więc zależność 1/N od w”*? jest prostą przecinającą oś odciętych w punkcie 
1/N, o nachyleniu, które wyznacza wartość k’ (wykres B na rycinie 13); 
4. zachodzą reakcje (107a) i (107b), k, = 0. Mamy wówczas 


1/N = (x + 1)/N,. (110) 


Efektywność zbierania nie zależy w tym przypadku od szybkości wirowania 
elektrody, tak więc wykres zależności 1/N od w”? jest linią prostą równoległą do 
osi rzędnych (wykres C na rycinie 13). Wartość parametru x zależy od potencjału 
elektrody dyskowej; 

5. W układzie zachodzą wszystkie rozpatrywane reakcje, w związku z czym 
równanie (108) nie może być uproszczone. Wykres zależności 1/N od w? jest 
prostą o nachyleniu zależnym od potencjału dysku (wykres D na rycinie 13). 

Mclntre (52) przedstawił kryteria pozwalające na odróżnienie kinetyki re- 
akcji równoległych (takich jak przedstawiono powyżej) od kinetyki reakcji ele- 
ktrodowych z katalizą heterogeniczną, której schemat przedstawiają poniższe 
równania: 


A >B (111a) 
M reakcje 

B 2C e (111b) 
ysku, 

pB“ A + qD (111o) 


gdzie reakcja (111c) jest powierzchniową reakcją rzędu p regeneracji składnika A, 
przy czym stała szybkości regeneracji wynosi k,. Powyższy schemat można prze- 
pisać tak, aby przedstawiał tok reakcji równoległej [52]: 


k 
AB (112a) 
BC (112b) 


A + p/p- Dlej> [a (p - D]D. (1120) 
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Ten schemat można interpretować, posługując się równaniem (108), przy czym 
parametr x definiowany jest równaniem 


gdzie 
n' = (p-8)/p, (114) 
oraz 
8=F:p:'k,' [c (O)]F/i,. (115) 


Na rycinie 14 przedstawiono wykres zależności 1/N od w? w przypadku, gdy 
w układzie zachodzi reakcja (111a) i (111c) dla p = 2 oraz przypadek, gdy k, = 0 
(linia przerywana). Dla reakcji katalitycznej x jest różne od zera w odróżnieniu 
od przypadku, gdy k, = = 0, stąd 1/N zależy od w”. Parametr x zależy od poten- 


: 00 01 02 03 04 05 0.6 0.7 08 09 10 


1t 


Ryc. 14. Schemat zależności 1/N od w' dla heterogenicznej katalitycznej reakcji elektrodowej II 

rzędu, w przypadku gdy w układzie zachodzą reakcje (111a) i (111c) (czyli dla potencjałów elektro- 

dy w zakresie I fali). AE=E-E?, gdzie E? jest potencjałem formalnym reakcji (111a). AE; > AE, 

>...AE;, linia przerywana oznacza przypadek, gdy stała szybkości heterogenicznej reakcji katali- 
tycznej (111c) wynosi zero 
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Ryc. 15. Schemat zależności 1/N od 1w!? dla heterogenicznej katalitycznej reakcji elektrodowej II 

rzędu w przypadku, gdy w układzie zachodzą reakcje (111a), (111b) i (1110) (czyli dla potencjałów 

elektrody w zakresie II fali). AF=E-E$, gdzie E$ jest potencjałem formalnym reakcji (111b). AE; 

>AE, >AE,. Linie przerywane oznaczają przypadek, gdy stała szybkości heterogenicznej reakcji 
katalitycznej (111c) wynosi zero 


cjału (113—115); stąd dla każdego potencjału inna krzywa zależności 1/N od w~”. 
Zakrzywienie powyższych wykresów odzwierciedla zależność parametru x od w. 
Na rycinie 15 przedstawiono z kolei wykresy zależności 1/N od w~” dla przypad- 
ku, gdy w układzie ząchodzą reakcje (111a, b ic) oraz p = 2. Linią przerywaną oz- 
naczono odpowiednie krzywe, gdy k, = 0 (czyli x = 0). Nachylenia krzywych 
w tych przypadkach rosną wraz ze wzrostem potencjału dysku. Dla ustalonego 
potencjału dysku wartości 1/N są większe w przypadku procesu katalitycznego, 
gdyż substancja B jest usuwana z układu poprzez elektrochemiczny proces 
(111b) oraz reakcję katalityczną (111c). Zakrzywienia na wykresach procesu kata- 
litycznego odzwierciedlają zmiany parametru n° wraz ze zmianą w. Przy potencja- 
łach bliskich potencjałowi prądu granicznego procesu B — C zachodzą relacje 
n'-> 1i x — 0, tak więc krzywe katalityczne będą dążyły do krzywych tego pro- 
cesu, w którym k, = 0. 
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5.2. WYZNACZANIE PARAMETRÓW KINETYCZNYCH REAKCJI ELEKTROCHEMICZNYCH 
METODĄ RRDE 


5.2.1. KRZYWB POLARYZACYJNE NA DYSKU I PIERŚCIENIU 


Analiza krzywych polaryzacyjnych umożliwia wyznaczenie potencjału półfali 
oraz liczby elektronów biorących udział w reakcji elektronowej na dysku. Analiza 
krzywych polaryzacyjnych na pierścieniu pozwala na ustalenie liczby związków 
rejestrowanych na pierścieniu oraz na wyznaczenie potencjałów półfali dla każde- 
go z tych związków. Krzywe woltamperometryczne na dysku i pierścieniu dostar- 
czają danych umożliwiających wyznaczenie wydajności prądowej procesu prze- 
biegającego na pierścieniu zgodnie z wyrażeniem 


Q, = |n,'ia,/N* n,*i,| * 100% . (116) 


5.2.2. WYZNACZANIE STAŁYCH SZYBKOŚCI REAKCJI ELEKTRODOWYCH METODĄ RRDE 
1. Proces elektrodowy dwustadyjny 
Rozpatrywany proces ma postać 
A+ne'eB+ne ŁC, (117a) 
przy czym produkt B jest rejestrowany na pierścieniu (ta reakcja przebiega szyb- 
ko): 
B +n,e >E. (117b) 
Wyrażenia opisujące zmiany stężenia substancji A i B na dysku mają następu- 
jącą postać: 


D,/B, > (c3 -c,(0)) = k, . (©). (188a) 


k,- c (0) =D, * cz(0)/65+ k, * ¢,(0). (118b) 


Z powyższych wyrażeń łatwo można otrzymać takie równania opisujące stęże- 
nia substancji A iB na powierzchni dysku: 


160 : Jarosław Penar 
c,(0) = cY(1 + k, -6,/D,); (119a) 


c,(0) = c, (0) - k/k,(1 + D, /k, * 8,)- (119b) 


Całkowity prąd dysku (i,) jest równy sumie prądów dwóch reakcji stadyjnych. 
Równanie opisujące ten prąd ma postać następującą: 


by = (ip) + (ip), =A F (n,*k, *c,(0)+ n, k, c,* cz(0)). (120) 
Podstawiając za cz(0) równanie (119b) mamy 


CA 


b= (m+n) (+D Aa) t nf(i + Dk, : 8,)], (121) 


gdzie (ip), jest granicznym dyfuzyjnym prądem dysku związanym z prze- 
pływem (n, + m) elektronów, tzn. (ip,), = A * (m, + m): F * Dy * c3/6A. Prąd ię re- 
jestrowany na pierścieniu stanowi część N, dyfuzyjno-konwekcyjnego prądu 
związanego z elektrodową generacją substancji B: 


ip = My *A *Ny*Dz *cz(0)/0, (122) 


Podstawiając za c,(0) równanie (119b) oraz (119a) otrzymamy następujące peł- 
ne wyrażenie opisujące prąd pierścienia: 


PY ng Ny (i wa j Dz: ôs 
n= GTW (1 +D,/E 3) (1+D,%, „) (123) 


Zatem 


_ Pat + r) 66) 


ini s (124) 


Podstawiając zaś za 5, równanie (15) otrzymamy 
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Mr (n, +n,) 1,61 k, D%.y% 4 
Zależność i,/i, od w”? jest liniowa, przy czym punkt przecięcia z osią odcię- 
tych wyznacza wartość n,, a z nachylenia można wyznaczyć stałą k,. Aby wyzna- 
czyć stałą k, i parametr n,, należy zależność (125) przedstawić w następującej 
formie: 


D,-4 D, 
EEE E +6, +, Be Hmm) eM ND» | (24a) 


= Wstawiając wyrażenie (15) za 6, i 5, , otrzymamy jawną zależność prawej stro- 

ny równania (124a) od w. Wykres powyższej zależności jest prostoliniowy, 
z wartością n, odciętą na osi [(ż,,), — ip]ię Oraz nachyleniem wprost proporcjonal- 
nym do k,. Tak więc, aby wyznaczyć wszystkie cztery parametry (k,, k,, n,, n,), na- 
leży przeprowadzić pomiary i, oraz i, dla co najmniej dwóch różnych szybkości 
wirowania elektrody. Wartość (i„,), może być wyznaczona bezpośrednio w eks- 
perymencie lub obliczona z odpowiednich danych. 

II. Proces elektrodowy równoległy 

Taki proces można przedstawić w następujący sposób: 


Mi, k 


y (126) 
= 


Wyrażenia opisujące zmiany natężenia substancji A i B na powierzchni elektro- 
dy mają postać 


(D,/8,) «(8 -c,(0)) = (k, +) ` c,(0), (127a) 
k, ` c,(0) = (D;/6,) * cz(0). (127b) 
Równanie opisujące natężenie prądu dysku ma następującą postać: 


i,=A*F(n,*k, +n, k) > c,(0). (128) 
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Z kolei równanie dotyczące natężenia prądu pierścienia ma następującą postać: 


ip = A+ F + tg: N,* Dy Cy (0)/8,,. (129) 


Laczac wyrazenia (127 a), (127 b), (128) oraz (129) otrzymamy 


+ pp _ ktm ką . 
infi = N me (130) 


oraz 

: —" n, DV w 
[GA — Via = nz No + Torn, Mg 3 Ny k, ik (131) 
gdzie (ip,), jest granicznym prądem dysku związanym z przeniesieniem n, ele- 
ktronów. Równania (130) i (131) umożliwiają — na podstawie pomiarów prądów 
dysku i pierścienia przy różnych szybkościach rotacji elektrody — na rozdzielenie 
reakcji (126) i wyznaczenie wartości k, i k, oddzielnie. 


5.3. WYZNACZANIE STAŁYCH SZYBKOŚCI REAKCJI HOMOGENICZNYCH METODĄ RRDE 


1. Reakcje homogeniczne I rzędu 
Rozpatrywany proces można przedstawić następująco: 


dysk A >B, 
roztwór B *» produkty, (132) 
pierścień B — mA, 


przy czym jeśli substancja B reaguje z nadmiarem innej substancji (np. C), to za- 
chodzi k, = k, © c%. Równanie różniczkowe opisujące rozkład stężenia substancji 
B na elektrodzie ma następującą postać [10]: 


D : (8%cy/02?) — v,* (dc,/dr) — v, * (dcB/d,) —k, 'cz=0. (133) 


To równanie można rozwiązać analitycznie tylko dla elektrod z wąskim pier- 
ścieniem i wąską przerwą izolacyjną [10]. Bierze się to stąd, że aby rozwiązać 
równanie (133) w strefie przerwy izolacyjnej i strefie pierścienia, należy zanie- 
dbać człon v, (dc,/dz) opisujący konwekcję normalną dla elektrody. Takie założe- 
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nie można przyjąć tylko dla elektrod o małych rozmiarach przerwy i pierścienia, 
gdyż wtedy decydujący wpływ na stężenie powierzchniowe substancji B w strefie 
przerwy i pierścienia ma konwekcja radialna. Jeśli mamy do czynienia z szeroką 
przerwą izolacyjną i/lub szerokim pierścieniem, to na stężenie powierzchniowe 
substancji B w tych strefach znaczny wpływ wywiera konwekcja w kierunku pro- 
stopadłym do elektrody (konwekcja normalna), tak więc człon trzeci w równaniu 
(133) nie może być wówczas zaniedbany. Przybliżone rozwiązanie równania 
(133) dla k,= 0 ma postać (przy zaniedbaniu konwekcji normalnej) 


N,”=-28 0 -28-8' +2: (0 + B) F/B’) + (B -3 2 TH - (BYE (a, (134) 
gdzie: 
a’ = 3: In(rJ/r,); P = 3 - ln (r./r,). 


Dla elektrody z wąskim pierścieniem i wąską przerwą. izolacyjną mamy 
a=(r/r,} —1—a oraz B = (ry/r,)* —1— B’, tak więc N, > N,”. Rozwiązaniem 
równania (133) dla k,=0 jest skomplikowane wyrażenie zawierające transformaty 
Leplace’a funkcji Airy. Analityczna retransformacja tego wyrażenia okazała się 
niemożliwa do przeprowadzenia [10], jednakże rozwinięcie funkcji Airy w wielo- 
miany, a następnie ich retransformacja pozwala na numeryczne rozwiązanie rów- 
nania (133). Dla parametru k < 1 (gdzie K= (k,/D)!? - (D/C), C= 0,51 - w”? - vt?) 
można otrzymać przybliżone, analityczne rozwiązanie równania (133). Ma ono 
następującą postać: 


NP=NP PYP: (1-A: v.)+0,5: Apts AZ ++ v.(P')8-2*4,'k2*T, (135) 


gdzie — dla x < 0,3 — v. = A; [1 — 0,372: (A, * K} + 0,146 - (A, ' k)*...], a dla K > 
>0,3 v. = K“ * tanh (A, ' x). Parametr v. jest bezwymiarowym stężeniem powie- 
rzchniowym wprowadzonym przez H ale’ a [53]. Funkcja 7, jest zależna od wy- 
miarów geometrycznych elektrody [10]. Ponieważ w praktyce rzadko spełniony 
jest warunek N,” =N; można, ze względu na trudności w konstrukcji elektrod z 
wąską przerwą i wąskim pierścieniem, wartość N,” „urealnić” poprzez odpowied- 
nie równanie korekcyjne [10] 


N, = (1+ (1-tanh(A,x)) : [(N/N”o”)-1)]:Ny” (136) 
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Powyższe wyrażenia wskazują, iż wartość N, zależy jedynie od wymiarów ele- 
ktrody oraz od wartości kinetycznego parametru «K. Stąd, wykreślając zależność N, 
= f(x) dla elektrody o danych wymiarach, można wyznaczać wartość k znając 
eksperymentalnie wyznaczone wartości N,. Gdy k > 1, wartość N, oblicza się z 
przybliżonego równania (137) [10] 


logN, — logN,= —a'x+ b — łogf1+ 1/(K+ C)]. (137) 


Współczynniki a, b i c zostały obliczone i stabelaryzowane w monografii A 1b e- 
ry'ego iHitchmana [10] dla częściej spotykanych w praktyce wartości r/r, i 
r,/r,. Przedstawiona powyżej metoda wyznaczania stałej szybkości reakcji homoge- 
nicznej I rzędu jest przydatna dla wartości parametru x z przedziału (0,2;5), co odpo- 
wiada wartościom stałej k, z przedziału (3-10*s', 10s") oraz wartościom k, (stałej re- 
akcji pseudopierwszego rzędu) z zakresu (3*102dm*/mol-s,107dm*/mol's). 

Powyższa metoda nie może być niestety stosowana w przypadku elektrody z 
szeroką przerwą izolacyjną i/lub szerokim pierścieniem. W takim przypadku le- 
piej jest stosować metodę wykorzystującą modelowanie numeryczne reakcji (132) 
przebiegającej na RRDE. Opis modelu oraz wyczerpujące omówienie założeń 
programu numerycznego wraz z metodyką pomiarów, a także porównanie tej me- 
tody z opisaną powyżej, można znaleźć w literaturze [54, 55, 56, 57, 58]. W sy- 
mulacji cyfrowej otrzymuje się krzywą roboczą zależności N, vs. XKT, gdzie 
XKT= «>. Z tej krzywej wyznacza się wartość parametru XKT dla wartości N, 
otrzymanych z pomiarów eksperymentalnych. Następnie wykreśla się zależność 
tak otrzymanych wartości XKT w funkcji c3/w. Nachylenie tak otrzymanej prostej 
jest wprost proporcjonalne do k,. 


II. Reakcje homogeniczne II rzędu 
W tym przypadku mamy do czynienia z następującymi procesami: 


dysk A + npe >B, 
roztwór B +C 2 produkty, (138) 
; pierścień B + nge > A. 


Równania dotyczące transportu masy mają następującą postać: 


Dp (0°c, /dz) — v, * (dc,/dz) — v, * (dc,/dr) — k, * cz: c.=0, (139a) 


Do: (©c J32) — v,- (c,/02) — v, (8c, /ar) — k, *cz*c.=0. © (139b) 
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Liniowo wzrastający prąd dysku generuje substancję B. Gdy reakcja homoge- 
niczna przebiega bardzo szybko (tzn. k, o), to stężenie substancji B na powie- 
rzchni dysku będzie równe zeru, dopóki prąd dysku nie będzie dostatecznie duży, 
aby wytworzyć taką ilość substancji B, która przewyższy ilość substancji C trans- 
portowanej do elektrody drogą konwekcji i dyfuzji. Wówczas 


M=A 'n,'F - D”? : C'? c0/1,288, (140) 
gdzie M to krytyczny prąd dysku. 


Gdy i,= M, to gęstość przepływu substancji B od elektrody jest równa gęstości 
przepływu substancji C do elektrody. Gdy i, > M, wówczas wytwarza się warstwa 
roztworu, w której zawartość substancji B znacznie przewyższa stężenie substan- 
cji C. Dalszy wzrost prądu powoduje rozszerzanie się warstwy z dominującą Za- 
wartością B w kierunku elektrody pierścieniowej. Gdy ta warstwa osiągnie we- 
wnętrzną krawędź pierścienia, to przez elektrodę zaczyna płynąć prąd i, (punkt a 
na rycinie 16). Wówczas 


i, = -M/[1-F(a)]. (141) 


Dalszy wzrost natężenia prądu dysku powoduje, że strefa dominacji substancji 
B wędruje wzdłuż pierścienia, aż cała powierzchnia pierścienia znajdzie się w tej 
strefie (punkt b na rycinie 16). Wówczas zmiany prądu dysku powodują liniowy 
przyrost prądu pierścienia według następującego równania: 


igy= N'i- M Pee (142) 


Wartość parametru M można zatem wyznaczyć przez ekstrapolację liniowej 
części krzywej ią, od i, do wartości ią, = 0. Rozwiązując równania (139a) i (139b) 
dla prądów dysku pomiędzy punktami a i b (Ryc. 16) — tzn. dla 0< B< B, gdzie B, 
= (r/r) — (ry/r,)?, a r, oznacza wartość r z przedziału (r, r,) — otrzymuje się war- 
tość efektywności zbierania elektrody N’, gdy strefa dominacji substancji B prze- 
suwa się wzdłuż pierścienia. Wówczas 


thira EJ 


N’ =1-F (wp) - Wang stn vena (143) 
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Prócz tego spełniona jest zależność 
Bu (1+a+ B)* 
Mii, = -PO E saapi lwia pagb bd (144) 


Powyższe dwa wyrażenia umożliwiają wyznaczanie wartości M bardziej pre- 
cyzyjne niż z krzywej ią, od i, [10,59]. Należy mianowicie z równania (143) wy- 
znaczyć wartość N” dla różnych B, w przedziale (0,8) oraz dla samych B, obliczyć 
wartość M/i, z równania (144). Następnie wykreśla się proste o nachyleniu N’ na 
krzywych eksperymentalnych ię; —i, (Ryc. 17). Punkty przecięcia prostych z krzy- 
wą eksperymentalną wyznaczają wartości (i,),,.„. Kąt nachylenia prostej (tp)... — 
(ip/M) jest równy M^. 

Gdy reakcja homogeniczna nie przebiega zbyt szybko, wtedy pomiędzy strefa- 
mi dominacji substancji B i C nie ma ostrej granicy. Część substancji B przechodzi 
do strefy dominacji substancji C i odwrotnie (Ryc. 18). To zjawisko może być wy- 


Ryc. 16. Krzywa miareczkowania w warstwie dyfuzyjnej na RRDE oraz rozkład stężeń substancji 
B i C pod elektrodą dla przebiegającej szybko reakcji homogenicznej II rzędu 
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Cig. dep 


(ip exp 
Ryc. 17. Schemat eksperymentalnej zależności natężenia prądu pierścienia od natężenia prądu dysku 
w przypadku reakcji homogenicznej II rzędu. Proste przecinające eksperymentalną krzywą igy vs ip 
obliczone zostały z równania (143) 


korzystane do pomiaru szybkości reakcji pomiędzy B i C, gdyż głębokość przeni- 
kania substancji B w obszar dominacji substancji C jest odwrotnie proporcjonalna 
do stałej szybkości tej reakcji. Głębokość penetracji łatwo jest mierzyć ustalając 
wartość prądu dysku w ten sposób, aby granica pomiędzy strefami B i C (tzn. ob- 
szar gdzie c= Cc znajduje się w r= r,, B= 0) znalazła się na wewnętrznej krawędzi 
pierścienia . Gdy reakcja homogeniczna przebiega bardzo szybko, przez elektrodę 
pierścieniową prąd nie płynie, ponieważ granica pomiędzy strefami B i C jest bar- 
dzo wąska (Ryc. 18). Jeśli jednak reakcja nie przebiega tak szybko, wówczas 
prąd pierścienia ma określoną wartość i to tym większą, im powolniejsza jest re- 
akcja. Zatem wielkość natężenia prądu pierścienia jest miarą głębokości przenika- 
nia substancji B w obszar dominacji substancji C. Gdy strefa reakcji znajduje się 
na wewnętrznej krawędzi pierścienia, to równanie opisujące natężenie prądu dys- 
ku — patrz równanie (144) dla B= 0 — ma następującą postać: 


li| =M/(1-F(a)). (141a) 


Natomiast tzw. kinetyczny prąd pierścienia (i,), — płynący przez elektrodę, gdy 
prąd dysku spełnia równanie (141a) — opisywany jest równaniem [60] 
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DYSK PIERŚCIEŃ 
Yb Y Yy 


Ryc. 18. Rozkład stężenia substancji B iC pod RRDE w czasie miareczkowania w warstwie dyfu- 
zyjnej w przypadku, gdy reakcja homogeniczna II rzędu nie przebiega szybko 


(i), = 0,21 + 2-72 -n-F + DO? >v i Az (145) 


Tak więc równanie dotyczące efektywności zbierania N’, — tzn. stosunek (i,),/ip 
— ma następującą postać: 


0,399 r3 D” [1—F(a)] 
N’,= RR a 0 OO > | (146) 


Z powyższych równań można wyznaczyć wartość k,. Analizując krzywe ię; od 
i, (jak na rycinie 17) należy wyznaczyć wartość M dla poszczególnych szybkości 
wirowania elektrody w w sposób opisany uprzednio, stosując zależności (143) i 
(144). Odpowiednie wartości M(w) wstawia się do równania (141a) otrzymując 
wartości (i,),(w) — kinetycznego prądu dysku. Z krzywej eksperymentalnej od- 
czytuje się wartości (i), dla danej wartości (ż,),(w), co daje w sumie wartości 
N*,(w). Wykreślając zależność N’, od w/c?, — zgodnie z równaniem (146) powinna 
być to linia prosta przechodząca przez początek układu współrzędnych — otrzy- 
muje się wartość k, z nachylenia prostej. Opisaną powyżej metodą można badać 
kinetykę reakcji o stałej szybkości k, do wartości 10°dm?/mol - s [61]. Trzeba jed- 
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nakże nadmienić, że zależność (146) w praktyce nie jest prostoliniowa. Dla w/c? 
—o występuje odchylenie od linii prostej, gdyż wartość N”, nie może dążyć do ©. 
Wartość k, wyznacza się więc wykreślając styczną do krzywej N’ («/c?. ) w zakre- 
sie niskich wartości œ/c% . Taka procedura obarczona jest zawsze dość sporym błę- 
dem. Prater iBard przedstawili metodę wyznaczania stałej szybkości k, opar- 
tą na technice symulacji cyfrowej [62]. Ta metoda polega na symulacyjnym wy- 
znaczeniu krzywych roboczych N, od bezwymiarowego parametru XKTC - m 
(gdzie (XKTC = k, : co: w+ v8 - DP -(0,51)%, m= ce/'c9) przy skoku potencjału 
na dysku do wartości prądu granicznego ż,,. Otrzymuje się w ten sposób (dla ele- 
ktrody o danych wymiarach) rodzinę krzywych; dla każdej wartoŚci parametru m 
jedna krzywa. Ustalając w eksperymencie wartość natężenia prądu dysku na war- 
tość graniczną i,,, należy rejestrować zmiany N,= t,,/ip,) przy zmianie wartości 
wic% Stałą szybkości k, otrzymuje się przez wyznaczenie z krzywych roboczych 
wartości XKTC - m (dla danego m) dla parametrów N, otrzymanych w pomiarach 


Ni 


RKTC*m 


Rys. 19. Symulacyjna krzywa robocza N, vs XKTC - m dla elektrody o wymiarach r,= 0,3 cm; r= 
=0,33 cm; r4= 0,42 cm, x — m= 0,001088, -m = 0,004354, -m= 0,006531, -m= 0,008710 
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i następnie na wykreśleniu funkcji XKTC * m od w/c?.. Nachylenie tak otrzymanej 
prostej jest wprost proporcjonalne do 1/k,. 

Inny sposób wyznaczania stałej k, (przy wykorzystaniu symulacji cyfrowej) 
podano w pracy [62]. W tej metodzie na elektrodę dyskową nakłada się prąd o 
wielkiej wartości, aby powierzchnia reakcji (Ryc. 18) znajdowała się na wewnę- 
trznej krawędzi pierścienia. Dla wartości N’, [61], otrzymanych z analizy do- 
świadczalnych krzywych ig vs. ip (Ryc. 17), znajduje się odpowiadające im warto- 
ści XKTC + m na roboczej krzywej symulacyjnej N, vs. XKTC - m. Dla danej ele- 
ktrody wartość N, jest taka sama dla układów z tą samą wartością XKTC : m, przy 
czym prąd dysku musi być tak dobrany, aby powierzchnia reakcji (Ryc. 18) znaj- 
dowała się w tej samej pozycji dla każdego z tych układów. Dzięki temu przy sta- 
łym prądzie, przy którym powierzchnia reakcji osiąga wewnętrzną krawędź pier- 
ścienia, otrzymuje się jedną krzywą roboczą (Ryc. 19) bez względu na wartość m. 
Wykreśla się następnie zależność XKTC - m od c?/w. Nachylenie prostej jest 
wprost proporcjonalne do 1/k,. Ta metoda ma kilka zalet w porównaniu z metodą 
Albery’ego i in. [61], opisaną powyżej. Charakteryzuje się większą do- 
kładnością i nie jest ograniczoną wielkością parametru c?/w oraz wymiarami ele- 
ktrody. Ma rownież tę zaletę, że obliczenia symulacyjne sprowadzają się do 
wygenerowania pojedynczej krzywej roboczej, co skraca czas obliczeń oraz upra- 
szcza procedurę analizy danych eksperymentalnych. Prócz tego dla wyznaczania 
wartości k, najkorzystniej jest „ustalić” pozycję powierzchni reakcji na wewnętrz- 
nej krawędzi pierścienia, gdyż pierścień „zbiera” wtedy jedynie tę część substan- 
cji B, która nie zdążyła przereagować z substancją C w obszarze jej dominacji. 
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SUMMARY 


The paper examines the usefulness of the rotating disc (RDE) and rotating ring-disc (RRDE) 
electrodes for investigations of the mechanism and of kinetics of electrode reactions. The study pre- 
sents in detail the description of the manner of determining the kinetic parameters of various types 
of electrode reactions. Mention is made of the digital simulation technique employed for the investi- 
gations with RRDE. 


